Absence of long-range superconducting correlations in the frustrated ^-filled band 
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We present many-body calculations of superconducting pair-pair correlations in the ground state 
of the half-filled band Hubbard model on large anisotropic triangular lattices. Our calculations cover 
nearly the complete range of anisotropies between the square and isotropic triangular lattice lim- 
its. We find that the superconducting pair-pair correlations decrease monotonically with increasing 
onsite Hubbard interaction U for inter-pair distances greater than nearest neighbor. For the large 
lattices of interest here the distance dependence of the correlations approaches that for noninter- 
acting electrons. Both these results are consistent with the absence of superconductivity in this 
model in the thermodynamic limit. We conclude that the effective i-filled band Hubbard model, 
suggested by many authors to be appropriate for the K-(BEDT-TTF)-based organic charge-transfer 
solids, does not explain the superconducting transition in these materials. 
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I. INTRODUCTION 



The two dimensional (2D) Hubbard model has been 
extensively investigated because at i-filling it can suc- 
cessfully describe the antiferromagnetic (AFM) phases 
found in many strongly-correlated materials. Since AFM 
phases often occur in materials displaying unconventional 
superconductivity (SC), such as the high-T c cuprates 
and the organic k-(BEDT-TTF) 2 X (hereafter k-(ET) 2 X) 
charge transfer solid (CTS) superconductors, it has fre- 
quently been suggested that some small modification of 
the model can yield a superconducting state where the 
residual AFM fluctuations mediate an attractive pairing 
interaction. In the case of the cuprates, this modification 
involves a change in the carrier concentration ( "doping" ); 
the doped 2D Hubbard model has been intensively inves- 
tigated with numerous analytic and numerical methods, 
but whether or not SC occurs within this model is still 
controversial. 



The AFM state in the 2D Hubbard model can also 
be destroyed at fixed carrier concentration by the in- 
troduction of lattice frustration. The model on the 
anisotropic triangular lattice (see Fig.[TJ has been used to 
describe the k-(ET) 2 X family of CTS, where SC occurs 
at fixed carrier density under application of moderate 
pressure. The ET layers here consist of strongly dimer- 
ized anisotropic triangular lattices, with the intradimer 
hopping integrals much larger than the interdimer ones. 
Each (ET) J dimer contains one hole carrier on the aver- 
age. This has been used to justify replacing each dimer 
unit cell with a single site, and the underlying ^-filled 
cation band with an effective i-filled bandi. 



We investigate ground state superconducting pair-pair 



correlations within the Hamiltonian, 

(ij),a [kl],a 

+ U y^nj^rii^. (1) 

i 

In Eq. [TJ c\ a creates an electron of spin a on site i and 

ni, a = c i(r Ci t(T . U is the on-site Hubbard interaction. We 
consider a square lattice with hopping integral t along 
x and y-directions and frustrating hopping t' along the 
x + y-direction (see Fig. [1]). The limits t'/t = and 
1 correspond to the square and the isotropic triangular 
lattices, respectively. All quantities with dimensions of 
energy will be expressed hereafter in units of t. We con- 
sider only the ^-filled band corresponding to an electron 
density per site p = 1. 

The non-superconducting phases of this model, shown 
schematically in Fig. [TJ are relatively well established. 
As t' is increased in strength, frustration destroys 
the q=(7r,7r) AFM ground state, replacing it with ei- 
ther a paramagnetic metallic (PM) state or a non- 
magnetic insulator (NMI) stated. Numerical calcula- 
tions on this model and the related model with two 
diagonal t' bonds in each plaquette have confirmed 
the presence of the PM, AFM, and NMI phaseaSriS. 
The NMI phase has been suggested as a candidate 
stated that explains the apparent quantum spin liq- 
uid (QSL) behavior seen in the strongly frustrated k- 
(ET) 2 Cu 2 (CN)3 1 . As the NMI phase has already been 
extensively investigate d 3 ' 4 ' 6 ' 9 ' 15 ' 16 , in the present work 
we will not consider the properties of this phase any fur- 
ther, but will rather focus on the possibility of SC within 
the model. 

Numerous mean-field theories have suggested that un- 
conventional SC occurs adjacent to AFM-PM phase 
boundarie o 7 ' 8 ' 20 " — . Similar superconducting states 
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FIG. 1: (color online) Lattice structure (see inset) and the 
ground state phase diagram of the | -filled band Hubbard 
model on the anisotropic triangular lattice for < t' < 1. 
Phases labeled are paramagnetic metallic (PM), Neel an- 
tiferromagnetic (AFM), and non-magnetic insulator (NMI). 
d x i_ y 2 superconductivity has been suggested to occur near 
the boundary of the AFM and PM phases (see text). Filled 
circles are finite-size scaled values for phase boundaries from 
PIRG calculations (see Section HV)| . The precise value of TJ c i 
at t' — 0.8 is known with less certainty. At t' = 0.8 the solid 
circle is an exact upper bound from the 4x4 lattice, while the 
dashed circle shows the expected finite-size scaled value. The 
phase boundary lines linking the points are only schematic 
guides to the eye. Dashed vertical lines indicate the parame- 
ter regions we investigate for superconductivity. 



have been suggested for the closely related Hubbard- 
Heisenberg model on the same latticed— . Because of 
the proximity of (7r,7r) AFM order, the suggested sym- 
metry of the SC order parameter is d x 2_ y 2. For t' ~ 1 
the magnetic ordering q shifts to (3-, 3-) corresponding 
to the 120° ordering found in the triangular lattice an- 
tiferromagnetic Heisenberg model, and consequently su- 
perconducting other order parameter symmetries have 
been suggested^ .. The estimated value of t' for the k-ET 
materials is however smaller than 1 (see below) 32 ' 33 , and 
also no evidence for 120° AFM order is found experimen- 
tally within the k-ET family 1 . 

Superconducting pair-pair correlations calculated with 
numerical methods going beyond mean field theory pro- 
vide a more accurate assessment of the presence of SC, 
provided finite-size effects can be adequately controlled. 
Two criteria must be satisfied to confirm SC within the 
model: (i) the superconducting pair-pair correlations 
must be enhanced over the U = values over at least 
a range of U, and (ii) the pair-pair correlations must ex- 
trapolate to a finite value at long inter-pair distances. 
We have previously calculated pair-pair correlations for 
the 4x4 lattice using exact diagonalizationi .. No en- 
hancement of the pair-pair correlations by U was found 
in these calculations, except for a trivial short-distance 
enhancement 1 ^ (see also below) . Our present work allows 
more careful analyses of the distance dependence of the 
pair-pair correlations, as well as the U-dependence of the 



longer-range components of these correlations, that were 
not possible within the earlier small cluster calculation. 

Pair-pair correlations for lattices comparable to 
those in the present work have also been calcu- 
lated using variational quantum Monte Carlo (VMC) 
methodsi 3 -^^^! 3 ^. VMC results however depend to 
a great degree on the choice of the variational wave- 
function, and there is considerable differences of opin- 
ion within the existing VMC literature. Clearly, studies 
of pair-pair correlations on large lattices, using many- 
body methods that do not depend on an a priori choice 
of the wavefunction are desirable. A candidate method 
for calculations of strongly-correlated systems is the re- 
cently developed Path Integral Renormalization Group 
(PIRG) method 3 ^. Like VMC, PIRG is also varia- 
tional and does not suffer from a fermion sign problem 
as do standard quantum Monte Carlo methods. Unlike 
VMC methods however, instead of an assumed functional 
form of the wavefunction, PIRG uses an unconstrained 
sum of Slater determinants that is optimized using a 
renormalization procedure 3 ^—. The NMI phase within 
Eq. Q]was first identified using PIRG 2,3 . Previous PIRG 
calculations 2 - - — investigated the metal-insulator tran- 
sition, AFM ordering, and properties of the NMI phase in 
detail, but did not discuss superconducting pair-pair cor- 
relations. Here we revisit the model with PIRG and cal- 
culate pair-pair correlations as a function of t' and U. As 
explained in Section Hill we use the most accurate version 
of the PIRG ground-state method, Quantum Projection- 
PIRG (QP-PIRG), which combines symmetries with the 
renormalization procedure^. As explained in Scction lHll 
we also performed an "annealing" procedure to help pre- 
vent the method from converging to local minima. 

While early tight-binding bandstructures calculated 
using the extended Hiickel method found some k-ET 
superconductors to have nearly isotropic triangular lat- 
tices with t' ~ 1, recent ab-initio methods have deter- 
mined that t' in the experimental systems lie within the 
range 0.4 < t' < 0.8 32 ' 33 . Importantly, in this range 
of anisotropy the 120° AFM order is not relevant. Fur- 
thermore, the AFM order is known experimentally to 
be of the conventional Neel patter n 40 ' 41 . Consequently, 
we limit our calculations to t 1 < 0.8. Specifically, we 
perform our calculations for three distinct t' = 0.2, 0.5, 
and 0.8, as shown in Fig. [TJ The two large t' values 
chosen bracket the estimated frustration in the k-ET 
superconductors 3 ^ 3 .. We choose a smaller t' = 0.2 in 
addition because it has been suggested in some studies 
that SC is present even in the weakly frustrated region of 
the phase diagra m 17 ' 26 ' 30 . The phase diagram (Fig. [I} is 
qualitatively different at t' — 0.2 because the NMI phase 
does not occur for t' < 0.52£&. While the estimate for 
degree of frustration is remarkably consistent between 
different ab-inito methods 3 ^ 3 ., the estimated value of U 
for k-(ET)2X is less certain. We therefore perform our 
calculations over a range of U starting from U = 0. 

The organization of the paper is as follows. In Section 
ITI1 we introduce definitions of the order parameters we 
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calculate. In Section Hill we describe the PIRG method. 
Section [IV] presents our data for t' = 0.2, 0.5, and 0.8, 
followed by discussions and conclusion in Section |V] 

II. ORDER PARAMETERS 

To determine whether SC is present near the metal- 
insulator (MI) transition, in addition to superconducting 
correlations we need order parameters to distinguish be- 
tween metallic and insulating phases. To locate the MI 
transition we use two different quantities. The first is the 
double occupancy D — (n^n^). As U increases, a dis- 
continuous decrease in D occurs at the MI transition^. In 
addition, we calculate the bond order Bij between sites 
% and j, 

% X C U*» • "'■ • ( 2 ) 

In the following we have labeled the bond order between 
sites linked by the t' bond as B' . This particular bond 
order is nonzero in the PM phase but tends to zero in 
the AFM Neel ordered phase because there electrons on 
sites connected by t' have parallel spin projections 1 ^. 

The operator A j = creates a singlet pair on lattice sites 
i and j: 

The pair-pair correlation function is defined as 

^(r) = j$>M<AtA, +rM >. (4) 

v 

In Eq. |4]the phase factor g(y) determines the symmetry 
of the superconducting order parameter. We consider 
two possible pairing symmetries, d x 2_ y 2 pairing (a = d in 
our nomenclature below) where g{v) = {+1,-1, +1,-1} 
and v{v) = {x,y, —x, —y}, and d xy pairing (a — xy) 
where g{y) = {+1, — 1, +1, —1} and r(v) — {x + y, —x + 
y,-x-y,x-y}. 

In the presence of superconducting long-range order, 
P a (r) for the ground state in the appropriate pairing 
channel must converge to a nonzero value for |r| — > 
oo. This is seen clearly for example in the 2D —U 
Hubbard modeler—. In the thermodynamic limit the 
long-distance limit of the pair-pair correlation function, 
P a (r oo) , is proportional 4 ^ to the square of the super- 
conducting order parameter, (A Q ) 2 oc \P( a (r — > oo)|. 
The magnitude of (A Q ) may further be used to set 
limits 4 ^ on the superconducting condensation energy, gap 
amplitude, and T c . 

III. METHOD 

The PIRG method has been previously used for 
a variety of strongly-correlated systems including the 
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FIG. 2: (color online) Comparison of PIRG and exact results 
for the 4x4 lattice with t' = 0.5. (a) Variance extrapolation 
of the ground state energy for [7 = 4. The symbol at AE=0 
is the exact energy; the line is a linear fit. (b) Pd(r) versus 
r. Open (filled) symbols are exact (PIRG) results. The inset 
shows relative percent error in Pd(r) as a function of r. (c) 
Pd{r) at the largest possible pair spacing on the finite lattice, 
Rd(fmax), as a function of U. Open (filled) symbols are exact 
(PIRG) results. The inset shows the variance extrapolation 
of Pd(r max ) for [7 = 4. 



2D Hubbard mode&, ^-filled frustrated 2D Hubbard 
models^&i&i^, and the ^-filled Hubbard model on the 
checkerboard lattic o 46 ' 47 . Details of the method are well 
described in these references. Here we discuss details 
of our PIRG implementation, and present comparisons 
with exact results which demonstrate the accuracy of the 
method for calculating pair-pair correlations. 

The PIRG method uses a basis of L Slater deter- 
minants, \4>i). For L — 1 this coincides with the 
Hartree-Fock (HF) approximation. In practical calcu- 
lations, maximum Us of a few hundred are used. The 
method is initialized with the L = 1 HF wavefunc- 
tion or a similar random starting wavefunction, and the 
PIRG renormalization procedure^ is used to minimize 
the energy by optimizing through the action of the oper- 
ator exp(— tH). One potential problem with the PIRG 
renormalization procedure is that the method can be- 
come trapped in a local minimum and not reach the true 
ground state 4 ^. Yoshioka et al. introduced a technique 
for avoiding local trapping by introducing global modi- 
fications to the wavefunction 4 ^. Following Ref. d(| we 
also introduced similar global modifications of the PIRG 
wavefunctions ("iteration A" in Ref. |46j), which mod- 
ify determinants in a global manner by acting on the 
wavefunction with exp[— tH] defined by a random set of 
Hubbard-Stratonovich variables. In addition, we also in- 
troduced updates to the \4>i) based on adding a random 
variation to the matrix elements of [4>i]j t k- The amplitude 
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of the variations is decreased systematically in a manner 
similar to simulated annealing. We found the addition 
of these two global updates to significantly improve the 
accuracy of the results. 

We also incorporated lattice and spin-parity symme- 
tries in the calculation^. Reference [38| introduced two 
different methods of using symmetry projection: (i) 
PIRG-QP, where symmetry projectors are applied to the 
ground state wavefunction after it has optimized using 
PIRG; and (ii) QP-PIRG, where symmetry projectors 
are applied at each step of the PIRG optimization. Here 
we have used the second more accurate of these two ap- 
proaches, QP-PIRG. The lattice symmetries we used in- 
cluded translation, inversion, and mirror-plane symme- 
tries, a total of 4JV symmetries where N is the number of 
lattice sites. We also applied the spin-parity projection 
operator after the PIRG process. An advantage of QP- 
PIRG is that much smaller basis sizes L can be used2&. 

Following reference [3?] we define the energy variance 
AE = ((H 2 ) - (H) 2 )/(H) 2 . AE is used to correct for 
the finite basis size L. For each set of parameters we 
performed the annealing and A iterations for successively 
larger basis sizes L. Each correlation function was then 
extrapolated to AE — by performing a linear fit. The 
error bars reported in our results are the standard errors 
estimated from the linear fit. The largest L used here for 
6x6 and 8x8 lattices was 256. The smallest L results we 
used in the fitting process depended on the value of U : 
for U < 2 we found that even L as small as 4 fit gave a 
good linear variance extrapolation, while for larger U we 
only used L > 16 results in the fit. 

In Fig. we compare results from our PIRG code with 
exact diagonalization results for the 4x4 lattice^. Here 
L of up to 256 were used. Fig. (2ja) shows the variance 
extrapolation of the energy for t' = 0.5. The extrapolated 
value for ground state energy is -15.037±0.002 compared 
to the exact ground state of -15.031. In Fig.^b) we plot 
the pair-pair correlation Pd{f) as a function of distance 
for £7=4. The inset shows the percent relative error in 
Pd{f) as a function of r. The maximum relative error is 
for r = r max =2\/2 and is smaller than 0.4%. Fig. [2][c) 
shows the d-wave correlation at the furthest distance, 
Pd{r mayi ) as a function of U for 4x4, t' = 0.5. The inset 
here shows the variance extrapolation for Pd(r max ) for 
U~A. Again, as in Fig. EJa) the extrapolation of the 
physical quantity is well within the statistical error. 

Our PIRG code was further verified against quan- 
tum Monte Carlo (QMC) results for larger lattices in 
the t' = limit where QMC does not suffer from the 
fermion sign problem at i-filling. For the 6x6 lat- 
tice, the QMC estimate for the ground state energy^ of 
Eq.[E]with U = 4 and t' = is E= -30.87±0.05. Previ- 
ous QP-PIRG calculations using lattice translations and 
spin-parity during the PIRG projection process, followed 
by a total-spin 5 = projection obtained^ E— -30.879. 
The extrapolated energy with our choice of AN lattice 
symmetries, spin-parity projection, and a maximum L of 
256 was almost identical, £'=-30.89±0.04. 





FIG. 3: (color online) PIRG results for t'=0.2. Squares and 
diamonds are for 6x6 and 8x8 lattices respectively, (a) dou- 
ble occupancy D (b) t bond orders (c) long distance d x 2_ y 2 
pair-pair correlation function, P<j(r*) as a function of U (see 
text). The inset shows Pd(r*) for the 8x8 lattice near the 
PM-AFM boundary, (d) Enhancement over non-interacting 
system, APd(r, U), for U = 2.75 as a function of distance. 
Circles here are exact 4x4 results. The inset shows the long- 
range part of APd(r, U) for the 8x8 lattice. In all panels lines 
are only guides to the eye. 



IV. RESULTS 

A. t' = 0.2 

At t' = 0.2 a single transition is expected between 
PM and AFM phases&Ifi. In Fig. E)[a)-(b) we plot D 
and B' as a function of U. The transition from the PM 
to an insulating phase is clearly seen as a discontinuous 
decrease in B' and D at U — U c . U c is only weakly 
size dependent at t' = 0.2 — for the 4x4 lattice^ JJ C = 
2.95±0.05, while for 6x6 and 8x8 lattices we found 2.75< 
U c <3.00. We estimate U c s» 2.7 in the thermodynamic 
limit. 

For all of the t' values we considered, we found that 
d x 2_ y 2 pair-pair correlations were of larger magnitude 
than d xy correlations (in Section [IV CI below we show an 
explicit comparison between the two). Fig. [3](c) shows 
the d x 2_ y 2 pair-pair correlations Pd{r*) as a function of 
U. The distance r* is defined as the next-to- furthest 
possible separation r between two lattice points on the 
finite lattice; r* = 2.24, 3.61, 5.00 for 4x4, 6x6, and 
8x8 lattices, respectively. Here we use r* rather than 
the furthest distance r max because of finite-size effects^ 
associated with r max . The 4x4 correlations are consid- 
erably larger in magnitude because of the larger r* on 
that lattice and we have not included them on Fig. [3]Jc). 
As seen in Fig. [3Jc), Pd(r*) has a tendency to decrease 
monotonically with U and is smaller at all nonzero U 
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FIG. 4: (color online) Short-distance (see text) d x i_ y i pair- 
pair correlation function, Pd{r — 0), as a function of U for 
t' = 0.2. Circles, squares, and diamonds are for 4 x 4, 6 x 6, 
and 8x8 lattices, respectively. Lines are only guides to the 
eye. 



compared to U 
ously. 



0. At U c Pd(r*) decreases discontinu- 



In addition to the [/-dependence, it is also impor- 
tant to examine the distance-dependence of pair-pair cor- 
relations. In Fig. [3Jd) we plot AP d (r,U), defined as 
AP d (r, U) = P d {r, U) - P d (r, U = 0), as a function of r 
for [7=2.75. Positive AP d {r, U) indicates enhanced pair- 
ing correlations over the noninteracting limit. We choose 
U = 2.75 in the PM state and close to the PM-AFM 
boundary where the greatest enhancement of pair-pair 
correlations from AFM fluctuations might be expected 
from prior work. Fig. [3]Jd) includes the exact 4x4 
AP d {r,U) as well. Our results in Figs. OUJc)-(d) show 
that as the system size increases, the long-range d x 2_ y 2 
pair-pair correlation function approaches that of nonin- 
teracting fermions. We have confirmed similar behavior 
of AP d (r, U), viz., absence of enhancement for other val- 
ues of U (not shown here) in either the PM or AFM 
regions. 

As seen in Fig.|3Jd), the only enhancement by U in the 
pairing correlations is at r — 0. The r — enhancement 
occurs because P d (r = 0) contains a component propor- 
tional to the nearest-neighbor spin-spin correlation func- 
tion; the enhancement of AFM order by U leads to an 
increase^ in P d (r = 0). In Fig. 2] we plot P d (r = 0) as 
a function of U. Precisely at U c there is a discontinuous 
increase in P d (r = 0), even as the system becomes semi- 
conducting, due to the increase in the magnitude of AFM 
spin-spin correlations. Importantly, only pair separations 
of r > 2 should be used to judge the enhancement of pair- 
ing correlations, because for r < 2 d x i_ y i pairs overlap 
on the lattice^. Here we find that P d (r) for distances 
beyond nearest-neighbor pair separation always decrease 
monotonically with increasing U. As we discuss further 
in Section [VJ the spurious increase of short-range corre- 
lations is the primary reason that mean-field calculations 
find SC near the MI transition. 
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FIG. 5: (color online) PIRG results for £'=0.5. Circles, 
squares and diamonds are for 4x4, 6x6 and 8x8 lattices re- 
spectively. Panels (a)-(c) are the same as Fig. EJa)-(c) except 
i'=0.5. (d) is the same as Fig. [3jd) except t'=0.5 and (7=4.5. 
In all panels lines are only guides to the eye. 



B. t' = 0.5 

Fig. E] shows D, B', and P d {r) for t' = 0.5. Not sur- 
prisingly, compared to t' = 0.2, U c \ depends here more 
strongly on lattice size, decreasing with increasing sys- 
tem size (U c i — 5.4 ± 0.1 and 4.6 ± 0.1 in the 6x6 and 
8x8 lattices, respectively). Previous PIRG calculations 
found U c i ~ 4.1 after performing finite-size scaling^. Our 
results are consistent with this value. 

In contradiction to t' — 0.2 (see Fig. d^b)), B' here 
is nonzero on the insulating side of the MI transition, 
suggesting that the nature of the insulating phase is dif- 
ferent. We have also calculated the spin structure factor 
S a {q) (not shown here). For U > U c i, a peak appears in 
S a {q) at q — (tt,tt). However, S(tt,tt)/N appears to ex- 
trapolate to zero as N —> oo, based on the three lattice 
sizes we have considered. This indicates that the sys- 
tem does not have long-range AFM order at t' = 0.5 for 
U > U c i, consistent with the NMI phase previously iden- 
tified in this parameter region^. Note that the larger B' 
in the NMI phase than in the AFM phase is also consis- 
tent with our previous exact diagonalization calculation 
(see Fig. 2 in Ref. [l(|) The properties of the NMI phase 
and the subsequent NMI-AFM transition at even larger 
U C 2 have both been extensively discussed in the literature 
before^— & 1 5 1 1 6 . Here therefore we focus on the strength 
of pair-pair correlations as a function of U and distance. 

Fig. 0[c) shows P d (r*) as a function of U. As at 
t' = 0.2, P d (r*) decreases monotonically with U. At 
U = U c i, Pd{r*) decreases discontinuously and is of very 
small magnitude in the NMI phase. The magnitude of 
P d (r) does not increase as U is increased further ap- 



6 




FIG. 6: (color online) PIRG results for t'=0.8. Circles, 
squares and diamonds are for 4x4, 6x6 and 8x8 lattices re- 
spectively. Panels (a)-(c) are the same as Fig. [5fa)-(c) except 
t =0.8. In (d) t' — 0.8 and U — 5.5. In all panels lines are 
only guides to the eye. 



proaching the AFM phase. Fig. [5fd) shows APd(r, U) 
as a function of distance for t'=0.5 for £7=4.5. As in 
Fig. IH at t' = 0.5 Pd(r = 0) increases at the MI transi- 
tion due to the increase in strength of nearest-neighbor 
AFM correlations, while long distance correlations are 
again weaker at nonzero U than at U = 0. 



C. t' = 0.8 

Our results for t' = 0.8, shown in Fig. [6j are similar 
to those for t' = 0.5. Here the PM region extends to 
somewhat larger 3 U. As at t' — 0.2 and t' — 0.5 there 
is no enhancement of the pairing correlations. Pd(r = 0) 
again shows an increase at the MI transition. From finite- 
size scaling we estimated the value for ?7 cl = 5.0 ± 0.3 
from our data. This value is identical to earlier results 3 . 

Several authors have suggested that the symmetry 
of the superconducting order parameter changes from 
d x 2_ y 2 to d xy or s + d xy in the region of the phase di- 
agram close to the isotropic triangular lattice limi t 25 ' 30 . 
In Fig. we compare the d x 2_ y 2 and d xy correlations 
for U — 4 as a function of r. Except at specific small 
r where pairs can overlap each other on the lattice^ -, 
we find that d x 2_ y 2 correlations are always stronger than 
d xy correlations. Plots of the P X y{ Y ) versus U also show 
a monotonic decrease with increasing U, and AP xy (r, U) 
similarly approaches zero for large r. 




r 



FIG. 7: (color online) Comparison of d x 2_ y 2 and d xy pair- 
pair correlations as a function of distance r at t' = 0.8 and 
U = 4. (a) 6x6 lattice (b) 8x8 lattice. In both panels, circles 
(squares) correspond to a = d (a — d xy ). Lines are only 
guides to the eye. 



V. DISCUSSION 

In Section iPvH we presented superconducting pair-pair 
correlations for the ground state of the ^-filled band Hub- 
bard model on the anisotropic triangular lattice calcu- 
lated using the PIRG method. Our main results are 
that (i) in all cases the superconducting pair-pair cor- 
relations at all finite U are clearly weaker than in the 
noninteracting limit, except for an enhancement of the 
very short-range component, and (ii) at large distances 
the distance dependence of the pair-pair correlations ap- 
proaches that of the noninteracting system. These re- 
sults, in conjunction with earlier exact diagonalization 
results 10 , which show exactly the same trends, strongly 
suggest that the superconductivity is not present in the 
model. Since many of the earlier works did find SC 
within the same model Hamiltonian, it is useful to com- 
pare these approaches and results with ours. Broadly 
speaking, two different kinds of methods had predicted 
SC within the triangular lattice Hubbard model, mean- 
field and the VMC. We discuss them separately. 
Mean-field approaches: In all cases, mean-field methods 
find a superconducting phase between the PM and AFM 
phases ^ 21 ' 22 . A NMI phase is found by some mean- 
field methods£&2i but not others 22 . However, there are 
further inconsistencies — for example, the paramagnetic 
insulating (PMI) phase in reference l2ll in some regions of 
the phase diagram unrealistically occurs at a temperature 
higher than that of the PM phase. In the NMI phase the 
nearest-neighbor AFM correlations are strong but AFM 
long-range order is not present. At fixed t', increasing U 
drives the system from the PM to NMI phases. Upon en- 
tering the NMI the d x 2_ y 2 correlation at r — increases 
in strength because the r = correlation is propor- 
tional to the nearest-neighbor spin-spin correlations. Our 
present results show that at the same time as this trivial 
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short-range correlation increases, the strength of longer- 
ranged correlations decrease greatly. Mean-field methods 
cannot capture these longer-ranged correlations and er- 
roneously extrapolate from the short-distance limit. 
VMC: Within VMC a functional form for the wavefunc- 
tion is assumed at the outset. The three different phases 
found for intermediate t' (PM, NMI, and AFM) require 
three different assumptions for the functional form of the 
wavefunction. This makes determining the true ground 
state behavior extremely difficult using VMC, especially 
near the phase boundaries. PIRG uses instead an un- 
constrained superposition of Slater determinants which 
does not suffer from these problems. The assumption of 
a functional form is a serious disadvantage as evidenced 
from the variety of results from different VMC studies 
which are not consistent with one another. For example, 
reference [13 did find SC in the model, while a later study 
by the same authors did not^. Liu et al. assume that the 
wavefunction is a Gutzwiller projected BCS function and 
find SC 2 ^. The appropriateness of such a wavefunction 
to describe SC is however a controversial assumption — 
Tocchio et al for example did not find SC within the same 
assumed wavefunction form 1 -. The occurrence of the 
NMI phase within VMC methods seems to be as much 
of a problem as within mean-field methods — Tocchio et 
al do not find evidence for NMI at i'=0.6 while PIRG 
does 3 -^. 

SC has also been found in several models that 
are closely related to the i-filled Hubbard model 
on the anisotropic triangular lattice. These include 
the Hubbard-Heisenberg mode l 17 ' 27 " —, and the Hub- 
bard model with two diagonal t' bonds per square 
plaquett o 26 ' 35 . While we cannot compare directly our 
PIRG results with these models, in nearly all cases the 
methods used to study these models are identical to those 
that have erroneously predicted SC within the present 
model. We have begun a reinvestigation of these models. 

We now return to the superconducting phase found in 
the organic CTS. Our results here cover the entire param- 
eter region (0.4 <t'< 0.8) thought to be appropriat e 32 ' 33 
for the K-phase CTS superconductors within the effective 
Hubbard model description for them, and clearly indi- 
cate that the ^-filled band Hubbard model is not suf- 
ficient to explain the occurrence of SC in /t-ET. It is 
important in this context to recall that in many CTS 
superconductors the insulator-superconductor transition 
is not from an AFM phase but from a different kind of 
exotic insulator. Examples include k-(ET) 2 Cu 2 (CN) 3 , 
which lacks long-range AFM order even at the low- 
est temperatures of measurement and has been consid- 
ered a QSL candidate 1 , but is superconducting under 
pressure^,, and other CTS superconductors where the in- 
sulating phase adjacent to superconductivity is nonmag- 
netic and charge-ordered (CO)^ 9 -. Once again mean-field 
theory, now within the -j-field band extended Hubbard 
Hamiltonian on anisotropic triangular lattices has sug- 
gested the possibility of a charge fluctuation mediated 
CO-SC transition^. Based on our present work, there 



are reasons to doubt mean-field approaches to SC within 
corrclated-electron models in general and these earlier re- 
sults should be checked through many-body calculations. 

The different natures of the unconventional semicon- 
ductors proximate to SC in the CTS confront theo- 
rists with a unique challenge. While mean-field the- 
ories suggest a different mechanism for each different 
semiconductor-superconductor transition, it appears un- 
likely to us that structurally similar materials, with iden- 
tical molecular components in some cases, should require 
different mechanisms for SC. Given how difficult a prob- 
lem correlated-electron SC has turned out to be we sug- 
gest that an alternate approach is to determine first how 
to construct a theoretical framework within which a uni- 
fied theory of SC begins to look feasible, and then to 
search for the same. We have recently shown how it 
may be possible to construct such a framework for the 
Q T g5i-54 In this picture, the k-(ET) 2 X and other dimer- 
ized CTS should be described in terms of the underlying 
2 -filled band as with the other CTS superconductor o 51 ' 52 . 
In the presence of strong dimerization and relatively 
weak frustration, AFM wins. Under increasing frustra- 
tion though, a transition occurs from AFM with uni- 
form charge density on each dimer to a charge-ordered 
paired singlet state we have termed the Paired Electron 
Crystal (PEC^i^. Experimental examples of the PEC 
in 2D CTS include /3-(meso-DMBEDT-TTF) 2 PF 6 (refer- 
ence [H and /3'-EtMe 3 P[Pd(dmit) 2 ] 2 (references MWk , 
which have precisely the same CO and bond patterns as 
in the PEC raodel 51 ^, and are superconducting under 
pressure. The application of pressure corresponds to a 
further increase in frustration and gives the possibility of 
a paired electron liquid superconductor—, a realization 
of the charged boson SC first proposed by Schafroth 59 . 
Although more work will be necessary to prove this, this 
theoretical approach has the advantage that it leads to a 
single model for correlated-electron SC in the CTS. Even 
more interestingly, we have pointed out that there exist 
several frustrated strongly correlated inorganic i-filled 
superconductors that can perhaps be described within 
the same model^— . 

Finally, the experimental observation of AFM^ in ex- 
panded fullerides A3C60 has led to the modeling of these 
compounds in terms of a 3D nondegenerate ^-filled band 
Hubbard model 61 . The threefold degeneracy of the low- 
est antibonding molecular orbitals in Cgo is removed by 
Jahn- Teller instability 6 -^. The observation of a spin- 
gap in the antiferromagnetic state validates the nonde- 
generate description 6 ^. The dynamic mean-field theory 
(DMFT) proposed for the AFM to SC transition in the 
fullerides within this 3D effective ^-filled band Hubbard 
model 6 -^ is however very similar to the DMFT theories 
of SC in the 2D CTS 7 -. Our results here suggest that a 
reexamination of the spin-fluctuation mechanism of SC 
in the fullerides may also be called for. 
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